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A simple theory of electromechanical transduction for single-charge-carrier double-layer electroac- 
tuators is developed, in which the ion distribution and curvature are mutually coupled. The obtained 
expressions for the dependence of curvature and charge accumulation on the applied voltage, as well 
as the electroactuation dynamics, are compared with literature data. The mechanical- or sensor- 
performance of such electroactuators appears to be determined by just three cumulative parameters, 
with all of their constituents measurable, permitting a scaling approach to their design. 



Electro-mechanical materials have long intrigued 
physicists. Such materials can serve as sensors and trans- 
ducers (where mechanical stimulation generates a volt- 
age) or can function as electroactuators where voltage 
generates a mechanical response (for review see [H-Q)- A 
class of such materials, called lonomeric Polymer Metal 
Composites (IPMC), is based on the motion of ions of one 
sign in an elastic polymer membrane enclosed between 
two bendable electrodes. Applied voltage causes mobile 
ions to accumulate near the oppositely charged electrode, 
forming an electrical double layer that expands the mem- 
brane and a depletion layer at the opposite electrode that 
contracts the membrane, causing bending (Fig. 1). 
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ation) for which it was straightforward to incorporate 
the stress-to-voltage feedback, as well as the volume of 
ions in electrical double layers. The analysis shows which 
physical characteristics determine the response and how 
they can be grouped into only two dimensionless cumu- 
lative parameters to describe the beam deflection under 
the applied voltage, with a third cumulative parameter 
determining the time scale of bending dynamics. This 
grouping reveals the interrelations between the physi- 
cal characteristics, instructs what should be measured to 
predict electroactuator performance, and makes possible 
a scaling approach to their design. 

Statics - The equilibrium properties of this system can 
be derived by minimising the free energy functional 
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FIG. 1. Electroactuators based on mobile cations: redistribu- 
tion of cations in the ionomer membrane between two flexible 
electrodes causes bending 0|. 

Early materials were swollen with polar solvents such 
as water, and were often single- ion conductors 0-01 ■ 
More recently, these volatile solvents were replaced with 
ionic liquids |7H10| to eliminate evaporation problems, 
but necessarily creating dual-ion conductors that often 
bend first forward and later backwards under a constant 
applied voltage 



11 



Sophisticated models for ionic ac- 
tuation with multiple sorts of charge carriers have been 
developed 12-1^. Those pioneering works, however, did 
not incorporate the feedback of mechanical deformation 
on the ionic charge distribution, which impedes their ap- 
plication for a description of the sensing function of elec- 
troactuators. The complexity of those models make them 
impractical for direct use in designing new actuators. 

We focus here on a simpler system, containing only 
one mobile charge carrier [l6| (for optimal electroactu- 
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where S and are the contact surface area and pro- 
jected area of the electrodes, respectively, fl7^; h is the 
half-separation of electrodes (positioned aX x — —h and 
X = h); e IS. the dielectric permittivity (assumed to be 
constant across the membrane); iIj{x) is the electrostatic 
mean-field potential distribution; p{x) — c+{x) — cq is 
the local excess number density of cations; K is the bulk 
modulus of the polymer matrix; v is the excess volume 
associated with each cation (including solvation shell, if 
present); k is the curvature of the actuator; E is the 
elastic modulus of the whole construction including the 
metallic plates, which in the first approximation can be 
considered as the volume weighted average of the Young's 
moduli of the metal and the polymer Il2ll : Cmax is the 
maximal local concentration of cations [l8j. 

The first term in Eq. [T]is the self-energy of the electric 
field plus that of the battery, keeping a defined potential 
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across the plates. The second term accounts for the elec- 
trostatic energy of all ions interacting with each other 
and with the external field. The third term (see E' 
is the entropy of the standard "lattice-gas model" [1 
including the volume of ions . The fourth term is the 
pressure- volume work that provides the coupling between 
mechanical bending and the ion concentration profile. In- 
deed, Kvp is the pressure exerted by ions and kx the 
elongation of the device due to curvature at a distance 
X away from the centre (a; > in extension; 2; < in 
compression). This effective pressure- volume work char- 
acterises ion-induced swelling. This term is proportional 
to the sum of first moments of the charge distributions at 
each electrode, reflecting the fact that the more polarized 
the charge distribution, the more favourable bending is 
and vice versa. The fifth term, outside the integral, is the 
energetic cost of bending an Euler-Bernoulli beam. 21 1 
Note that the beam is assumed to be much longer than 
its width, precluding Gaussian curvature with only 1- 
dimensional bending. 

All in all, Eqs. Il|2l describe several competing physical 
effects: the electrostatic and entropic costs of building up 
ions near one electrode and depleting them at the other. 
Bending relieves the energy of charge accumulation but 
the extent of bending is inhibited by the mechanical en- 
ergy needed to deform the material. 

Minimising the free energy functional (Eqs. Ill2p with 
respect to p, constraining the total number of ions, gives 
the electrochemical potential of the mobile ions /i+ = 
eip + kBThi[c+/{cma,x — c+)] — KvKx. By equating to 
its value in the bulk, 22] /iq = fcsrin[co/(c,nax~co)jj. we 
obtain the cation distribution (Fermi-like if k = 0) 11 
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Consequently, minimisation of Eqs. [T] and [2] gives a (di- 
mensionless) modified Poisson-Fermi equation 
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with dimensionless variables X = x/Id, y = eip/ksT — 
PkXId where p = Kv/ksT, ion crowding parameter 
7 = co/cmax, Debye length = (Itt/bCq) and Bjer- 
rum length Is = e^/eksT. Equation |4] is solved with 
boundary conditions y' = y = Q&tx = Q, correspond- 
ing to vanishing electric field and potential in the bulk. 
These conditions are true for widely separated electrodes, 
so that the electrical double layer does not overlap the 
depletion layer, maintaining electroneutrality in the bulk; 
warranted unless the membrane is nanoscale thin. One 
can then combine two semi-infinite solutions correspond- 
ing to the anode and cathode. Taking the first integral 
of Eq. 131 gives the expression for the electric field: 
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Minimising Eq. [T] with respect to k with H = h/lo 
and p = (c+/co) — 1, gives the dimensionless curvature: 
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The approximate equality is justified as the range of the 
charge density variation is much shorter compared to the 
separation between the electrodes 23|. Hence, x ^ Q 
where Q is the charge per unit cross-section area stored 
in both the double layer and depletion layer, consistent 
with experimental results [s']. The integrals in Eq. |5]can 
then be evaluated using the Gauss law ,24] which relates 
surface charge densities at each of the electrodes to elec- 
tric inductions. Since the electrical induction (Eq. [S]) 
contains itself x through the definition of y, a transcen- 
dental equation for x emerges: 
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Here a, the dimensionless electroactuation number, is the 
key cumulative parameter for the actuation strength [25| . 
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The response vanishes when the membrane is soft or the 
plates are rigid {K/E — ^ 0). Here, the dimensionless po- 
tential drops V, Vi and V2 are, respectively, the total 
voltage, and the potential drops near the positive and 
negative electrodes, all normalized to ksT/e, obtained 
via equating the induction, thus the surface charge den- 
sity, on both electrodes (this transforms the system into 
a constant voltage ensemble): 
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Expanding Eq. [7] in the low voltage regime gives Vi 
V2 w V/2 and x w a-^/ (1 — 7)^. In the opposite large 
voltage limit, Vi ~ {1 — "f)V and x w a^S (1 — 7) V. 
Figure 3 shows the full numerical solution of Eq. [71 where 
these two limits are recovered. This solution is replicated 
by the interpolation formula 
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matching the two limits. As seen from Eq. [TUl the 
crossover from linear to non-linear regimes occurs at 
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FIG. 2. Bending reduces the charge separation between the 
electrodes. Charge density profiles (Eq. [3]) are plotted as 
functions of the distance from the corresponding electrode 
within the proposed self-consistent theory at applied 4 Volts 
and for a = 1.3 x 10"^ 7 = 0.25 and S/So = 5. 2H = 2hl-^ , 
is the membrane thickness measured in units of Debye length. 



V ^ 8. At 300K, this corresponds to 0.2 Volts (see 
Fig. 3). This crossover has been observed experimen- 



tally in single-mobile-charge systems [2 
rationalized in ref 28[. Rearranging Eq 
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a simple law for the treatment of experimental data: 
(V/x)'^ = (1 — 7) (1 + which gives a straight line 

plotting (V"/x)^ vs V. Experimental tests of this law and 
its 1/8 ratio between the slope and intercept are needed. 

To illustrate how the full solution of Eq. [7] works, con- 
sider an example: 1 fim thick Nafion of volume fraction 
0.6, with volume fraction of mobile cationic counterions 
0.1 and ethylene carbonate solvent 0.3 with "typical" 
physical characteristics (dielectric constant e ~ 10 [29| . 
ion excess volume v — 0.3 nm'^ [s^l, ion crowding param- 
eter 7 = 0.25 313 and K/E = 13 for Nafion 

Figure [2] shows that the effect of curvature on the con- 
centration profile is to reduce ion polarization at the elec- 
trodes. This is expected, as bending relieves the steric 
and electrostatic strain induced by the accumulation and 
depletion of mobile ions near the electrodes. Bending 
opens extra volume close to the negative electrode, low- 
ering the local concentration of cations. Conversely at 
the positive electrode, bending decreases the extent of 
the depletion layer. 

Figure 3 demonstrates that larger electrode/membrane 
contact area and larger volume of the charge carrier en- 
hance bending. More charge is stored for higher sur- 
face area electrodes and crowding in the double layer 
increases with ion size [35|. To achieve a substantial 
curvature of k = Icm"^ at an applied voltage of 4V, 
assuming 7 — 0.25, a/h must be greater than 0.03cm~^. 
Using the above physical characteristics, the separation 
between the electrodes must be less than lO^m. 

This simple model has no hysteresis or dissipation ef- 
fects; the electroactuator response is fully reversible. Ex- 
perimentally, hysteresis may occur due to non-idealities 
such as Joule heating, leakage of solvent outside the elec- 
trode, irreversible mechanical deformation of the ionic 




FIG. 3. Dimensionless curvature x as a function of ap- 
plied voltage. The main panel demonstrates the effect of 
ion crowding parameter 7 at fixed a. The insets show two 
of the many ways to vary a; the effects of roughness fac- 
tor 5/5*0 and of the effective volume of ion for 7 = 0.25, 
5/ So ~ 1 and v = ^[BMimi+ • The volumes of ions used are 
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hydrated ion volumes Vi^i+ = 0.23 nm"', Vc^+ = 0.15 nm'' [34|]. 
Other physical characteristics are the same as Figure 2. 

polymer metal composite fsol, Is?! and the current from 
trace water hydrolysis. Nanoporous electrodes [sl, 
are also not considered; surface roughness is assumed to 
occur on a scale much larger than the Debye length. 

Dynamics - The simplest results can be obtained as- 
suming: (i) the mechanical response time of the polymer 
is much faster than the charging of the double layer, (ii) 
the height of the electrode roughness is much smaller 
than the electrode-electrode separation and (iii) the di- 
mensions of all channels are much larger than both the 
hydrated ion size and the Debye length. In this case, 
the response dynamics are solely determined by the mi- 
gration of ions through the bulk and an equilibrium-like 
charging of the electrical double layers. The theory of 
double layer charging in single-mobile-charge-carrier sys- 
tems has been considered in detail 1^ 4^ . Applying this 
to the case of electro-actuation, in response to a dimen- 
sionless step potential Uo{= eV^kBT) imposed at i = 0, 
we find an equation for >c{t) I4ll|, 
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Here r = RqCo] Rq = 2h/{Soa) = 2hkBT/{SoCoe^D) is 
the resistance of the bulk (cr is the ionic conductivity and 
D is the diffusion coefficient of the ion); Co — eS/{8nlj:)) 
is the linear Debye capacitance, giving 
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Equation 11 is solved analytically in the limits of linear 
response (C/q < 8) and strongly nonlinear response: 
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FIG. 4. The main panel shows fitting of experimental data 
taken from refs [3l with the non-linear response case of 
Eq. 13. There, a rectangular Nation membrane of dimension 
4x 35 X 0.2mm was coated with a platinum electrode and then 
treated with a thin gold overlayer. The counterion was Li"'" 
with glycerol as solvent. A IV applied potential (C/o — 40) 
induced bending. Assuming 7 = 0.25, the two-parameter fit 
to Eq. 13 yields a = 5.6 x 10~^ and r = 28s. The two insets 
show the effects of the electrode contact surface area and ion 
size on actuation response, from a 4V step voltage (see text 
and the Figure 3 caption for other parameters used). 



The two forms coincide at short times, where the initial 
rate of bending aUo/r is independent of 7, since crowding 
of ions requires time to develop. The maximum final 
bending (Eq. [TOl) scales as a^/T^^ which is maximal at 
Co = Cinax/2 (or 7 = 1/2). Hence, a vital prediction of 
our model is that single-ion actuators could be improved 
by roughly doubling the ion exchange capacity |42i of 
Nafion ionomers. The relation between excess ion volume 
and actuation rate can be obtained by noting that in the 
linear response regime, >f ~ vt/Ro with Rq v^^^ [43| . 
This gives the counter-intuitive prediction >c ^ v^/^t that 
the actuation becomes more rapid with larger ions! In the 
main part of Figure 4, the full x{t) curve, calculated via 



Eq. 13, is compared with experimental data 

Conclusion - A self-consistent theory of single-ion con- 
ducting electroactuators is presented, in which the curva- 
ture is considered on the same footing as the concentra- 
tion profiles, modelling both actuation and sensing. The 
model assumes a strong coupling of ion build-up and de- 
pletion with local volume that forces bending. The the- 
ory reveals the roles of all physical characteristics but 
most importantly, they group themselves into three cu- 
mulative parameters: (i) dimensionless response coeffi- 
cient - the electroactuation number a (Eq. [8]), (ii) ion 
crowding parameter 7 and (iii) the relaxation time r for 
dynamics (Eq. 12). These three parameters have clear 
constituents, each of which can be independently mea- 
sured or estimated. Other results are as follows: 

- At large voltages, both curvature and charge accumu- 
lation increase as owing to ion crowding. 

- For a given voltage, increasing the surface roughness, 
volume of ion and decreasing the electrode-electrode sep- 



aration all increase bending. 

- Plotting experimental data as {V/x)^ vs. V is pre- 
dicted (Eq. [TU|) to give a straight line with a slope 8 
times smaller than the intercept. 

- The dynamic response is described by simple analytical 
formulae (Eqs. 12 and 13). 

- Electroactuation number a allows a scaling approach 
for rational design of electroactuators. [Jl 
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